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The ground-state properties of the spin- 1/2 antiferromagnetic Heisenberg model with spatially 
anisotropic couplings on a square lattice are investigated by a spin-rotation-invariant Green's- 
function approach and by Lanczos diagonalizations on lattices up to 36 sites supplemented by 
finite-size scaling. We focus on the anisotropy-driven transition from the Neel state to a para- 
magnetic state with antiferromagnetic short-range order and on the spatial dependence of spin 
correlation functions. Our principal result is that a rather sharp crossover in the magnetic behavior 
occurs at the coupUng ratio -Ro ^ 0.2 {R = Jy/ Jx)- 
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Low-dimensional spin systems described by quantum 
antiferromagnetic (AFM) Heisenberg models have at- 
tracted increasing attention, in particular with respect 
to the crossover from one to two dimensions [1-5]. In 
this context, one basic problem is the influence of spa- 
tial anisotropy on the staggered magnetization m in the 
ground state of the two-dimensional (2D) spin-1/2 AFM 
Heisenberg model 
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Here R — Jy/ Jx (throughout we set Jx = 1), and {i,j)x.y 
denote nearest neighbors (NN) along the x-, y-directions 
on a square lattice. In the model (1) there occurs a 
transition from a long-range ordered (LRO) Neel anti- 
fcrromagnct at i? = 1 (m = 0.3074 [6]) to a spin liq- 
uid with pronounced AFM short-range order (SRO) at 
i? = (m = 0). The essential question is whether the 
critical coupling ratio Rc for the order-disorder transi- 
tion is zero or has a finite value. Whereas the ordi- 
nary spin-wave theory [1] gives Rc = 0.034 and the 
one-loop renormalization-group analysis of an effective 
spatially anisotropic nonlinear sigma model results in 
Rc ~ 0.047 [7], the RPA spin-wave theories [8,9] and 
the chain mean-field approaches [1,10] yield Rc — 0. 
From the Fade approximants to Ising expansions for m, 
Rc ^ 0.02 was suggested [3]. As stated by Affleck et 
al. [3], renormalization-group arguments do not give a 
definite answer whether or not the system orders for arbi- 
trarily weak interchain coupling Jy. Recently, Sandvik [5] 
developed a multi-chain mean-field theory and presented 
quantum Monte Carlo (QMC) results for an effective 3- 
chain model and for full 2D lattices (at i? = 0.1 and 
R — 1) which provide strong evidence for Rc — 0. 

Experimentally, in the quasi-lD antiferromagnets 
Sr2Cu03 and Ca2Cu03 a very small ordered moment 
and Neel temperature Tpf were observed and found to 
decrease with increasing anisotropy [11]. Based on a de- 
tailed estimate of the exchange integrals, Rosner et al. [9] 
obtained ratios of the order of i? ~ 10~^ and concluded 
that the previous theories overestimate both m and Tjv. 



Motivated by this unsettled situation, in this Letter we 
examine the order-disorder transition and the spin cor- 
relation functions of arbitrary range in the 2D spatially 
anisotropic Heisenberg model by an analytical approach 
based on the Green's-function projection method and by 
a finite-size scaling analysis of exact diagonalization (ED) 
data. Our results indicate a rather sharp crossover in the 
spatial dependence of the spin correlation functions in the 
model (1) at the coupling ratio Rq ~ 0.2. To this end, we 
extend the spin-rotation-invariant theory of Refs. [12,13], 
which yields a good description of spin correlation func- 
tions in the isotropic model (i? = 1), to the anisotropic 
case. 

To determine the dynamic spin susceptibility 
(q,^^) = —{{S^; SZq}}uj by the projection method, 
we consider the two-time retarded matrix Green's func- 
tion ((A;At))ij in a generalized mean-field approxima- 
tion [13] 



((A;At)), = [. 
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with the moments M ^ ([A, At]) and M' ^ {[lA, At]). 

We choose the two-operator basis A = {S^,iS^)^ and 
get 



x"^ (q, ^) = 



M, 



(1) 
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where M^^^ ^ {[iS+ , S^]) is given by 

M^i) = -4Ci,o(l - cosfc) - 4i?Co,i(l - cosg^,) (4) 

with Cr = (SqS^) = Cn^rrn and r = + mey. The 
Neel ordering in the model (1) is reflected by the closure 
of the spectrum gap at Q = (tt, tt) and by the condensa- 
tion of that mode. Thus, at T = 0, we have 



N ^ 2u„ 



(5) 
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The condensation part C results from (5) with r = 
employing the sum rule Co,o = \- Then the staggered 
magnetization is calculated as 
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l^(SoSr)e-^«- = ^C. 



(6) 



To obtain the spectrum lu^ in thc approximation 
WqiSq , we take the site-representation and decouple the 

products of three spin operators in —S^' along NN se- 
quences, introducing vertex parameters in thc spirit 
of the decoupling scheme proposed by Shimahara and 
Takada [12]: 



al'y{S+Sr)St+a,{StSr)S+ 
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Here, af and a\ are attached to NN correlation functions 
along the x- and y-directions, respectively, and a2 is as- 
sociated with longer ranged correlation functions. We 

obtain 



-[1 + 2a^^Cuo + 4Ra\Co,i 
+2a2(C2,o + 2RCui)] cosq^ 

-R[R + AalCifl + 2i?af Co,i 
+2a2(i?Co,2 + 2Ci,i)] cofiqy 

+2alCifl cos 2qx + 2R^afCo,i cos 2qy 

+AR{a^Ci^o + a\Co.i) cosqx cosqy . 



(8) 



Note that our scheme preserves the rotational symmetry 
in spin space, i.e., x (^1,^^) = x(q,w) = ^x'^-{q,u}). 

Considering thc uniform static spin susceptibility x = 
limq^o -^q^V^i^qj thc ratio of thc anisotropic fmictions 



M^^) and < 
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Iq^ must be isotropic in the limit 



0. That is, the condition 



(9) 



has to be fulfilled. The theory has nine quantities to 
be determined self-consistently (five correlation functions 
in Wq, m, and three vertex parameters) and eight self- 
consistency equations (six Eqs. (5) including Co,o = 1/2, 
thc LRO condition ujq = 0, and Eq. (9)). If there is 
no LRO, we have luq > 0, and the number of quanti- 
ties and equations is reduced by one. As an additional 
condition wc treat three cases for the choice of the free 
a parameter: (i) In case A, we adjust thc ground-state 
energy per site e{R) taken from the Ising-expansion re- 
sults by Affleck et al. [3] (Fig.l). (ii) In case B, we fix the 
uniform susceptibility x(^) by a linear interpolation be- 
tween the exact value x(0) = I/tt^ and the QMC result 
X(l) = 0.0446 [14] (star and open circle in the inset of 
Fig.l, respectively), which is justified by the ED data of 
Ref. [2]. (iii) In case C, we fit the free parameter to thc 
QMC results for m(i?) of the 3-chain mean-field theory 
by Sandvik [5] (Fig.2). 

Figure 1 compares thc ground-state energy s{R) ob- 
tained by different approaches. Here the ED data 
is extrapolated using the scaling law s{N) — e{oo) oc 



j^-{R+2)/2 ^ with an effective exponent agreeing with the 
known values in the ID [15] and 2D [16] cases (note that 
the 2D scaling exponent —3/2 seems to be more appropri- 
ate down to i? ~ 0.25). In order to reduce cluster-shape 
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FIG. 1.: i?-dependence of the ground-state energy per 
site e and the uniform static spin susceptibility x (inset). 

effects, we use two different groups of clusters of equal 
symmetry. Our ED data agrees with the Ising-expansion 
results by Affleck et al. [3] and, in particular, with the 
Bethe value e(0) = -0.4431. In the region R < 0.25, thc 
ground-state energy of case B approximately reproduces 
the exact data. Equivalently, at i? < 0.25 the suscepti- 
bilities x(^) in the cases A and B nearly coincide (see 
inset of Fig.l). 

Our results for the order parameter m plotted in Fig.2 
indicate an order-disorder transition at the critical ratio 
i?c — 0.24 (note that the cases A and B give very similar 
results for R < 0.25, cf. also Fig.l). The linear decrease 
of m with R in the region 0.6 < R < 1 in case B is as- 
cribed to the crude (linear) interpolation of x{R) in that 
case. To analyze the ED data for m?{N) calculated by 
Eq. (6), we use the finite-size scaling arguments given 
in Refs. [16] and [17] for the 2D (frustrated) Heisenberg 
antiferromagnet and fit the data to the scaling relation 
w?{N) — m^(oo) cx 7V~^/^ (inset of Fig.2). Thereby, as 
noted in Ref. [17], the extrapolated (A^ — > oo) values for 
m slightly depend on the factor in front of the sum in 
Eq. (6): A^"^ vs. [N + 2)-'^. The finite-size extrapola- 
tion of TO for A = 16, 36 depicted in Fig.l (1/A prefac- 
tor) and calculated with an 1/(A -h 2) prefactor agree, 
at i? = 1, with Ref. [17]. We find a transition to a spin- 
liquid phase at Rc ~ 0.24 or Rc ~ 0.18 depending on 
the chosen prefactor. Let us emphasize the coincidence 
of the critical coupling ratios obtained in cases A and B, 
and by Lanczos diagonalizations, where thc critical ra- 
tio Rc turns out to be nearly identical for both cluster 
sequences. To illustrate the finite-size scaling in more de- 
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FIG. 2.: Staggered magnetization m for different values 
of R. The short-daslied lines are fits to the ED data 

using A'' = 16, 36 (squares) and tilted N = 8, 16, 32 
(diamonds) clusters with periodic boundary conditions, 
respectively. The inset demonstrates the deviation of the 
ED results from the iV~^/^ scaling law at i? < 0.3 (sohd 
lines are least-squares fits). 

tail, in the inset of Fig.2 we show the least-squares fits of 
our ED results. The N~'^/'^ scaling [16], when applied to 
the small systems (up to 6 x 6), breaks down for R < 0.25. 
As recently pointed out by Sandvik [5,18], at i? < 0.2 the 
QMC simulations reveal an anomalous (non-monotonic) 
scaling behavior for square lattices, where very large sys- 
tems are needed to access the N~^/'^ law indicative of 
Neel order. 

To capture the behavior of m{R) with R^ = 0, let us 
consider case C. The results for e{R) and x{R) shown in 
Fig.l strongly deviate, in the physically most interesting 
region of R, from those obtained in cases A and B. In 
particular, at R 0.2, £{R) and x(-^) exhibit a strange 
maximum. 

Our findings indicate a rather sharp change in the mag- 
netic properties at i?o — 0.2. Note that Affleck et al. [3] 
found a poor convergence of their Fade approximants just 
below R < 0.2. To shed more light into the crossover at 
Ro, we calculate the dependence on R of the two-spin 
correlation functions depicted in Fig. 3. The NN corre- 
lation functions in cases A and B reasonably agree with 
the ED data for a 6 x 6 system, whereas Co,i in case C 
notably deviates from those results for i? < 0.25. That is, 
in the R region, where a finite magnetization is incorpo- 
rated (case C) as compared to the m = solution in the 
cases A and B, the SRO is described less well. The same 



FIG. 3.: NN and longer ranged (inset) spin-correlation 
functions versus R. Theoretical results are compared 
with ED data obtained for different cluster sizes/shapes. 
In the inset, at i? = 0, stars denote the extrapolated val- 
ues of the ED data for chains up to 32 sites using the 
scaling relation Cn,o{N) — C„,o(oo) oc N~'^. 

qualitative behavior is found for the correlation functions 
Ci,i and C2.0 (see inset of Fig. 3): At i? < 0.25, the results 
in cases A and B arc; in better agreement with ED data 
than the curve in case C. Note that the ratio Co,i/Ci,o 
determining the anisotropy in the spin-wave velocities ac- 
cording to Eq. (9) does not show a decoupling transition 
(Co,i/Ci,o = 0) at a finite R value, contrary to the sug- 
gestion by Parola et al. [2] . 

From our results we suggest the following characteris- 
tics of the crossover at i?o — 0.2. In our self-consistent 
Green's-function theory the short ranged magnetic cor- 
relations, in particular the fluctuations of the correlation 
functions about their LRO limit, are described rather 
well. On the other hand, the suppression of LRO by 
quantum spin fluctuations is overestimated. Moreover, 
if the correlation functions decrease slowly (subexponen- 
tially) towards their long-range limit, our theory does 
not capture this behavior. Equally, the ED data for the 
small systems cannot be described by the LRO scaling 
law N~'^/'^. Our results indicate that this happens for 
R < Ro — 0.2. Contrary, at R > Rq the decrease of the 
correlation functions with distance is sufficiently strong 
so that the LRO limit is nearly reached at relatively short 
distances. This yields a reasonable scaling behavior of 
our ED data (inset of Fig.2). To sum up, our results 
indicate that the crossover at i?o — 0.2 may be accompa- 
nied by a rather sharp change in the spatial dependence 
of spin correlation functions from strong to weak decrease 
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with distance. At R < Rq, the system starts to sense the 
ID behavior, characterized by an algebraic decrease of 
the correlation functions. 

Finally, in Fig. 4 we plot the spin- wave spectrum in the 
ordered state for the cases A (i? > i?o — 0.24) and C 
{R > 0) scaled by = c^(i?. = l)/V2, where = 
1.36 and Z^' = 1.58. For i? = 1, in both cases we have 
= 4^2(1 - 7^) with 7q = {cosqx + cosqy)/2. With 
decreasing i?, the anisotropy in Wq develops gradually, 
where in case C the ID limit Ug^ = 2y/—afCifi\ sin^'a,! is 
approached. However, compared with the exact Bethe- 
ansatz result LOg^ = ^ | sin \ , the spin- wave energy at 
<?a; = tt/S is too high by a factor of about three. 

To conclude, the strong evidence for the sharp 
crossover at the interchain coupling ratio Rq ~ 0.2 found 
by analytical and numerical methods should be corrob- 
orated by further QMC studies, especially on the dis- 
tance dependence of spin correlation functions. This is 
of particular interest in developing appropriate theories 
to explain the magnetism in quasi-lD quantum spin sys- 
tems [9,11]. 
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FIG. 4.: Spin- wave dispersion cjq of the 2D spatially 
anisotropic AFM Heisenberg model along the major sym- 
metry directions of the Brillouin zone. 
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